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Abstract 

In this paper the amorphous/solid to disorder liquid structural phase transitions 
of an anomalous confined fluid is analyzed using their local fractal dimension. 
The model is a system of particles interacting through a two length scales po¬ 
tentials confined by two infinite plates. In the bulk, this fluid exhibit water-like 
anomalies and under confinement forms layers of particles. The particle dis¬ 
tributions of them, present different arrangements related to amorphous/solid 
phases. Here only the contact layer is analyzed through fractal singularity spec¬ 
trum. At high densities the structural transition its quantify by the order de¬ 
gree to determine the phases affected by the confinement. This mapping shows 
that the system as the temperature increased, the fractal dimension decreases, 
which is consistent with the behavior studying in such systems. This result sug¬ 
gests that under thermodynamic perturbations, an anomalous confined liquid, 
presents different phase transitions achieving be characterized by its fractality. 
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1. Introduction 


The characterization of the phases present in complex system is not trivial. 
Usually it depends in identifying the correct order parameter of the structure. 
For instance, the structural transformation by thermal or mechanical perturba¬ 
tion of fibrous, dendritic, or colloidal configurations, formed by aggregation or 
reaction processes have been quantified from different measures of complexity. 
One of these measures is the fractal dimension. The fractality is a geometrical, 
topological and structural property present in many natural, physics or simu¬ 
lated complex systems [l], Q, Q] • In many cases a fractal structure results from 


the kinetic aggregation of a group of particles or from the reaction processes 


HQ, 13 1 


between them [1|, |4j, 15J1 in a process that resembles a very slow nucleation and 
growth of mechanism. This is visually manifested by different final distributions 
of the particles that entails universal properties, and also influence the physical, 
biological or chemical properties on the system Q, Q, Q, 7, Q. Q . These distri¬ 
butions can be also quantified by the mass fractal dimension that is a measure 


Q[n[Q. 


employed to quantify the different structural phase transition 

In a number of systems a single fractal dimension is unable to capture the 
full complexity of the system. The multi-fractal spectrum describes the scaling 
correlations, coexisting in the dynamical evolution of the system, at different 
length scales of observation. It is employed to provide a description of the aggre¬ 
gation kinetics. It also gives the information about how the new phase reaches 
the equilibrium state H ibl [uil . In this way the measure of the microscopic 
multi-scale structure through the local fractal dimension, is an important tool 
to identify the macroscopic state of the system, which influence the physical 
emergent properties. 

Then, the fractal dimensions can be used as an additional tool for charac¬ 
terizing the complex phases emerging from phase transitions. This strategy was 


n . 




employed in the study of rheological fluids [15|,[l6(, granular materials 7|, 17, 18] 


. 


n 


magnetic wall domains in boracite [19| and other complex systems [111, 20 


2M- 


In the case of rheological systems, the final structures obtained by magnetic 
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particles dispersed in mineral oils and perturbed by magnetic fields show dif- 


H- 


ferent degrees of order that were quantified by its mass fractal dimension 
This result was also checked experimentally [l5|, |l7|. This structural transfor¬ 


mations can be analyzed by glass transition approach [23, [23j, [24] and is possible 


to identify the liquid-glass and the liquid-crystal phase transition. 

All the examples cited above in which the fractal analysis were used to 
identify new phases were complex systems. Would this framework also be useful 
for describing phases in simple systems? In a simple fluid the thermodynamic 
and dynamic behavior is governed by the molecular length scale. This is the 
case of the rare gases, diatomic and triatomic molecules. In a complex fluid, 
the thermodynamic and dynamic properties are governed not by the atomistic 
length scales but by a mesoscopic scale that arises from the competition of 
the multi-scale molecular forces. These systems include colloidal suspensions, 
gels and polymer blends. Due to the complexity of the competition forces, 
complex fluids can be considered homogeneous at the macroscopic scale, but 
are disordered at the microscopic scale, and possess structure at an intermediate 
scale. This is the reason why the multi-fractal spectrum employed to analyze 
the structural transformation by mechanical perturbations can be applied in 
those complex fluids as well. 

Water, even though a very simple triatomic molecule, is not a simple liquid. 
It is an anomalous material showing a number of thermodynamic and dynamic 
anomalies [25]. The most familiar anomaly is its increasing density with temper¬ 


ature, at ambient pressure, up to 4°C. Above this temperature water behaves 
as a normal liquid and density decreases as temperature rises. Experiments for 
water allow to locate the line of temperatures of maximum density (TMD) in 
the pressure-temperature plane. Below TMD, density decreases with decreasing 
temperature, differently from the behavior of the majority of fluids, for which 
density increases on lowering temperature [261. 


In addition to the thermodynamic and dynamic anomalies, water exhibits 
many solid phases. Several coexistence lines separate the multiple solid phases. 
Thus, the energy landscape associated to the crystalline phases presents a num- 
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ber of sharp valleys with very low energies. The temperature and pressure 
ranges at which each one of these sharp valleys displays lowest energy values 
define the stable phase in that region of the phase diagram. Those valleys of the 
energy landscape that never achieve the lowest energy correspond to the amor¬ 
phous configurations. When confined within plates, the fluid energy landscape 
becomes even more complex. The anomalous fluid forms layers and the_system 


shows a transition from three layers to a two layers structure 


27 


28j. 29(- Us¬ 


ing nanotubes, the same transition appears and it is associated with a dynamic 
transition from a normal to superflow regime fl [dll . [32I . [dill . At low tempera¬ 
tures and high degree of confinement the contact layers melt, while the central 
layer stays liquid. The contact layer form a variety of liquid, liquid crystal and 
solid structures [34 [. 

Recently a model for describing the anomalous behavior of water were stud¬ 


ied under confinement 


27]. These studies indicate that the confined system 


exhibits at the wall two dimensional phases not present in the bulk system 


M 


35 
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36l 37]. While the existence of the phases is identified clearly by the 


instabilities of the density versus pressure phase diagram, the nature of the new 
structures, tested with the radial distribution function ]27| and with the trans¬ 


lational order parameter 


28 
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35 


30 


37j, it is still unclear. Particularly 


the system presents phases that change continuously to very different struc¬ 
tures without phase transition while other phases change through a first order 
transition. These two scenarios can not be identified by the translational order 
parameter analysis and need further understanding. 

In this work we explore the idea that the fractal analysis can provide infor¬ 
mation of the structure and phase behavior of anomalous fluids, like water. In 
this context we study the phase behavior of a water-like model confined within 
plates. The pressure versus temperature phase diagram, of this fluid is analyzed 
in the framework of the multi-fractal spectrum and within this framework the 
different phases are identified. 

The paper is organized as follows: in the section II the model is introduced; 
in the section III the methods are presented; the results are given in the section 
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IV and our final conclusions are presented in the section V. 


2. The Model 


The water-like fluid is composed by N spherical particles of effective diameter 
a that interact through a core-softened potential of two length scales, namely 


U{n 3 ) 


= 4 


12 


aexp 


1 


rij ~ r 0 


(i) 


where rij = |r, — fj | is the distance between two fluid particles i and j. The first 
term is a standard 12-6 Lennard-Jones (LJ) potential [3] and the second one 
is a Gaussian well centered at ro, with depth a and width co- The parameters 
used in this work are a = 5.0, Co = 1.0 and ro/cr = 0.7, that result in a 
potential with two length scales, one around r l:j = n « 1.2<r and the other at 
Tij = V 2 ~ 2cr [39]. The potential is shown in the Figure Oja). This model does 


not exhibit the directionality or explicitly hydrogen bonds as present in water, 
however it captures the competition between open and close water tetramers 
through the competition between the two length scales. As a result, the pressure 
versus temperature phase diagram shows a re gion where density, diffusion and 

3 43] and in confined 


structural properties are anomalous in bulk 


systems [27], [32 


33, 
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In our system the particles are confined between plates in ^-direction. All 
simulations were done with plates of size L x = L y = L = 20<r, separated by a 
fixed distance L z . A schematic depiction of the system is given in the Fig. |T] 
(b). The interaction between the fluid particles and the molecules in the plates 
is purely repulsive and it is given by the Weeks-Chandler-Andersen (WCA) j3] 
potential, namely 


U 


WCA 


0 ij) = 


4e 

[M 12 - 

f~) 6 

- 

\ z a ) 

\ Zi * J 


Zij < 2 1/,0 CT 


Zij > 2 1//6 cr . 


( 2 ) 


where is the distance between the plates at j position and the ^-coordinate 
of the fluid particle i. 
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Figure 1: (a) Two length scale interaction potential and (b) schematic depiction of the confined 
system. The particles are confined between two flat and smooth plates, separated by a distance 
L z . The energy and length scales are given by e and cr, respectively. 


3. Methods 


3.1. The simulation details 


The physical quantities are shown in reduced units 



r 


* 


r ^ * 3 

— and p = pa , 
a 


(3) 
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for distance and density of particles, respectively, and 


Pm = 


Pi a 3 


and T* = 


k B T 


( 4 ) 


for the pressure and temperature, respectively. The symbol * will be omitted 
in order to simplify the discussion of the paper. 

The molecular dynamic simulations were performed at NVT-constant. The 
systems have 500 particles confined between two parallel flat plates in ^-direction. 
The plates have thickness of a, area of L 2 and are separated by a fixed distance 
L z . The value of L was 20 in all the simulations and the density of the system 
was changed varying the value of L z , from 4.3 to lO.Oer. 

The repulsive interaction with the plates creates an exclusion region in the z 
direction, consequently the total density is corrected to an effective density [49, 
[ 5 ^ . Then the distance between the plates can be approach by L ze ps L z — a, 
resulting in an effective density of p = N/(L ze L 2 ). In x and y directions periodic 
boundary conditions were employed. 

The temperature of the systems was fixed using the Nose-Hoover heat- 
bath 5lJ,l52| with a coupling parameter Q = 2. Simulations were performed for 


the temperatures ranging from T = 0.050 to T = 0.450. This choice of temper¬ 
atures were made on basis of the pressure versus temperature phase diagram of 
the bulk system. The initial configuration was generated placing the fluid par¬ 
ticles randomly between the plates. The equations of motion were integrated 
using the velocity Verlet algorithm, with time step of St = 0.001. The systems 
were equilibrated with 5 x 10 5 steps followed by 1 x 10 s steps for the production 
of the thermodynamic averages. A particle-particle interaction cutoff radius of 
r cut/o' = 3.5 was used. 


3.2. The multifractal spectrum 

We used a standard box counting method, also called the capacity of the 


set 


3 £ 


M 


551 ], to calculate the fractal dimension on the two-dimension con¬ 


figurational structure of each final equilibrium state, on the layer of the fluid 
confined. The images are sets of two-dimensional final stages of the simulation 
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process. Each one present a particle configuration for fix densities and diverse 
values of temperature. The procedure to measure the fractal dimension is the 
follow: the image is binarized by a high contrast treatment, leaving the particles 
black and the space between them white. A grid of four random positions cover 
the entire image with a decreasing size of e as the length of the box. The scaling 
law to relate the number of particles and the size of the boxes follow the relation 

N ~ £~ d « (5) 


where e acquired successively smaller values of length until the minimum value 
of £q and N(e) are the number of cubes required to coverall the set. The fractal 
dimension by the box counting method is given by 

In N (e) 


D q = lim 

e->0ln(£ 0 /£) 


( 6 ) 


To describe all the statistical properties by the local fractal dimension 
we used the generalized box counting method { 57 !. I 58 . 


59 


53 


561 , 


601 defined as 


1 lnl(q,£) 

D q = -hm ——— 

1 — q e-s -0 ln{£ 0 /£) 


where 


N(e) 

= Et^W 5 

i—1 


(7) 


( 8 ) 


1 . H 


as 


P? 


l,q 


We used the scaling exponent defined by Halsey et al. 
where a can take a width range of values measuring different regions of the set. 
The spectrum generated by an infinite set of dimension D q = Dq, D\, D 2 , ... 
measure the scaling structure as a function of the local pattern density. If 
q=0 the generalized fractal dimension represent the classic fractal dimension, 
it means that Df = D q=0 [61]. As the image is divided into pieces of size £, 
it suggested that the number of times that a in Pi q takes a value between a' 
and dcJ defined as do!p[o!)£~P a ^ where f(a') is a continuous function. As q 
represents different scaling indices, we can define 

N(e) r 

I(q, e) = E = j da'p(a')£-^ + ^ (9) 


2=1 
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cq is the Lipschitz-Holder exponent, that characterizes the singularity strength 
in the ith box. The factor cq allows to quantify the distribution of complexity 
in an spatial location. The multifractal is a set of overlapping self-similar con¬ 
figurations. In that way, we used the scaling relationship taking into account 
f(a) as a function to cover a length scales of observations. Defining the number 
of boxes as a function of the Lipschitz-Holder exponent N(a ), can be related to 
the box size e as 

N(a) ~ £- /(a) (10) 


The multi-fractal spectrum show a line of consecutive points for Q > 0 that 
start on the left side of the spectrum climbing up to the maximum value. Then 
the values for Q < 0, represented in the spectrum for a dotted line, on the right 
side begins to descend. The maximum value corresponds to Q = 0, which is 
equal to the box counting dimension. To obtain the multi-fractal spectrum we 
use the plugin FracLac for ImageJ[62|. Basically D q is the variation of mass 


as a function of £ in the image, and give us the behavior as a power series of 
e sizes distorting them by an exponent q. We select the case of Df = D q= o as 
the parameter of order in the images. In the plugin we select four grid positions 
that cover the total image, and the mode scaled series was selected to see the 
singularity spectrum results. The final configuration of the simulations, present 
a particle distribution of particles in black. The parameters of the program were 
calibrate for this kind of images. 


4. Results 


The temperature versus density pha se diagram of this system was obtained 


27] showing the various two dimensional 


by molecular dynamic simulations 
phases present at the contact layer. 

The existence of phase transitions in the analysis of Bordin et al. 


271 was 


obtained by computing the density versus pressure isochores and observing in¬ 
stabilities related to first-order phase transitions. The general classification 


employed for by Bordin et al 


27] for defining the different phases of this system 
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took into account: the radial distribution function, the particles mobility and a 
snapshot of the system. It is interesting to observe that large structural changes 
can occur without phase transition. In order to understand how very different 
structures are connected here we explore the fractal dimension of each one of 
these structures. 

Now let us explore each density region of the phase diagram. Figure [2] (a) 
illustrates the temperature versus total density phase diagram showing the dif¬ 


ferent phases present at the fluid contact layer 


271. Figure [5] (b) show selected 


snapshots of the different phases. At this high density region at low temper¬ 
atures a liquid-crystal phase (I) coexists with a solid hexagonal phase (II). At 
higher values of the temperature both phases become liquid and increasingly 
disordered and for p c 3 = 0.321 and T c3 = 0.45 the coexistence disappears at 
critical point. It is important to notice that no phase transition is observed be¬ 
tween the liquid-crystal phase: phase I, and the liquid phases: phases III and V. 
Also no phase transition is observed between the solid hexagonal phase, phase 
II, and the other two liquid phases, phases IV and VI. 

Here we complement this analysis by identifying the different phases with 
the correlation between each final equilibrium state of the phase transition and 
the degree of ordering. 

The starting point for our analysis are the snapshots of the system. Fig. [2] (b) 
shows the pictures for the temperatures T = 0.075 (points I and II), T = 0.200 
(points III and IV) and T = 0.350 (points V and VI). These figures represent 
the final equilibrium state in which the specific structure exists. To make more 
evident this evolution, we calculated the fractal dimension of each layer at those 
different densities and temperatures. This order parameter gives a quantitative 
classification for these phase transitions. 

Figure EKa) presents the maximum values obtained in the singularity spec¬ 
trum as a function of total density, p 1 while figure [3|b) shows it as a function of 
the separation between plates, L z , for the three values of the temperature. In 
all the cases the fractal dimension increases with the increasing of confinement. 
Since the fractal dimension can be used as measure of the degree of order l.ll5|. 
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Figure 2: (a) Temperature versus total density phase diagram for high densities. The star 
represents the critical point at T c = 0.450 and p c = 0.321. The solid lines delimitate the 
coexistence region between liquid-crystal and hexagonal solid phases. The triangles, squares 
and circles are isochores at T = 0.075, 0.200 and 0.350, respectively, and densities pijuv ~ 
0.260, pu,iv,Vi ~ 0.370. The corresponding snapshots of the contact layer are shown in (b). 
More details in the text. 



Figure 3: (a) Fractal dimension as function of (a) total density and (b) separation of plates 
for transition at high densities. 

our results confirm that for each temperatures the order increases with density. 
It is interesting to notice that, the slopes of the constant temperature lines are 
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almost the same for all the temperatures analyzed, indicating that even though 
the structures change with temperature, the difference in order between the low 
density and high density structures, change proportionally. In addition, the 
values of the fractal dimension of the phases I, III and V are very similar, while 
the values for phases II, IV and VI are almost the same. This could explain 
why the transformation between them requires no phase transition. 



II 



IV 


i&WUi 



Figure 4: (a) Temperature versus total density phase diagram for intermediate densities. The 
star represents the critical point at T c = 0.20 and p c = 0.212. The solid lines delimitate 
the coexistence region between dimeric liquid and liquid-crystal phases. The snapshots of the 
contact layer are shown in (b) for the points I and II at T = 0.075 and for points III and IV 
at T = 0.100, with densities piju ~ 0.185 and pujv ~ 0.228. 


Fig. IDJa) illustrates the temperature versus density phase diagram of the 
liquid at the contact layer of the confined system in a region of intermediate 
densities 27]. At low temperatures a liquid-crystal phase (II) coexists with a 


structured liquid phase (I) made of dimers. When the temperature is raised both 
the liquid-crystal and the structured dimeric liquid becomes disordered, phases 
(IV) and (III) in the Fig.HJb) respectively. As the temperature is increased even 
further the coexistence ends at a critical point at T c 2 = 0.20 and p c 2 = 0.212. It 
is important to notice that even though quite different in structure it is possible 
to go from the phase I to phase III and from phase II to phase IV with no phase 
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transition. 

The snapshots of the final states of the system at the contact layer for this 
intermediate range of densities are depicted in the Figure |UJb). The figure 
illustrates the states I and II at T = 0.075 and for III and IV T = 0.100, 
with densities prjrr ~ 0.185 and pujv ~ 0.228, respectively. Employing these 



Figure 5: (a) Fractal dimension as function of (a) total density and (b) separation of plates 
for transition at intermediate densities. 


snapshots, the fractal dimensions of these configurations were computed. The 
Figure [5] shows the degree of ordering evolution as a function of total density and 
of the distance L z for this intermediate region of densities. Similarly to what 
happens at high densities shown in Fig. [3] the increase in the confinement and 
the decrease of the distance between the wall leads to an increase in the order. 
The phases I and III show almost the same value for the fractal dimension. 
The same occurs for the phases II and IV. This result support the idea that 
the transformation from one phase to another structurally very different with 
no phase transition might be related to our findings that they share the same 
fractal structure. 

Figure [SJa) shows the temperature versus density phase diagram of the liq¬ 


uid at the contact layer of the confined system in a region of low densities 


27 [ 


showing liquid configurations at coexistence. The coexistence between the two 
liquid phases ends at a critical point at T c i = 0.125 and p c l = 0.158. The snap- 
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Figure 6: (a) Temperature versus total density phase diagram for low densities. The star 
represents the critical point at T c = 0.125 and p c = 0.158. The solid lines delimitate the 
coexistence region between disordered liquid and dimeric liquid phases. The snapshots of the 
contact layer are shown in (b) for the points I and II at T = 0.075, with densities pj « 0.149 
and pn ~ 0.169. 



Figure 7: Fractal dimension as function of (a) density for T = 0.075 and (b) separation 
between plates. 


shots of representative phases at the contact layer are shown in the Figure [6jb). 
Phases I and II, with T = 0.075 and densities pi ~ 0.149 and pn ss 0.169 are 
illustrated, respectively. 

These equilibrium configurations were employed to compute the fractal di¬ 
mension of the system. In the Fig. [TJa), we show the fractal dimension as 
function of density, while in the Figure 0(b) the fractal dimension is computed 
as a function of the separation between plates. At this temperature, the disorder 
degree varies from Df = 1.6152 (p rss 0.150) to Df = 1.6511 (p « 0.170). Even 
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all the configurations being liquid states, we identify different degree of disorder 
in each case. 



Figure 8: Order evolution as function of L z for all the cases studied. 


In all the densities analyzed the fractal analysis confirms the results obtained 
with molecular dynamic simulations. Now in order to check if the different phase 
transitions show drastic differences in the evolution of the fractal dimension, the 
degree of ordering is compared. Figure [8] illustrates the fractal dimension as a 
function of the distance between the plates, L~ for all the cases studied here. 
The fractal dimension not only increases with the density but also a drastic 
change of slope is observed at very high densities. This region of pressures and 
temperatures represents the solid to liquid phase transitions where a high degree 
of ordering is expected confirming the suggestions of the molecular dynamic 
simulations. 

5. Conclusion 

In this work we have explored the use of the fractal dimension to analyze 
the phases present in an anomalous fluid confined by repulsive walls. 

The degree of configurational order-disorder of the confined liquid was an¬ 
alyzed using the fractal spectrum approach through image analysis. We found 
that different phases separated by phase transitions show a very different frac¬ 
tal dimension that increases with the increasing order of the structure of the 
phase. Complementary to this result we also found that phases that are not 
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separated by phase transitions show a very similar fractal dimension. This re¬ 
sult shade some light in the odd possibility of a continuous transition between 
two structurally very different phases. 
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